Current-driven domain wall motion with spin Hall effect: Reduction of threshold 

current density 
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We theoretically study the current-driven domain wall motion in the presence of both the spin 
Hall effect and an extrinsic pinning potential. The spin Hall effect mainly affects the damping ratio 
of the domain wall precession in the pinning potential. When the pinning potential is not too strong, 
this results in a significant reduction of a threshold current density for the depinning of a domain 
wall with certain polarity. We also propose one way to distinguish the spin Hall effect induced 
spin-transfer torque from the one induced by the Rashba spin-orbit coupling experimentally. 
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An electric control of magnetization by the spin- 
transfer torque (STT) Q is under intense investiga- 
tion due to its device application potentials. 0, 0] An 
essential prerequisite for the STT is to generate a spin 
flow of electrons. A simple intuitive method is to uti- 
lize spin-polarizcd conduction electrons in a ferromag- 
netic layer (FM). Recently, a very different way 
to generate a spin flow was proposed, |9l-[l2l| which uti- 
lizes a pure spin current generated by the spin Hall effect 
(SHE). When a charge current flows in a nonmag- 
netic layer (NM) with strong atomic spin-orbit coupling, 
the SHE generates a pure spin current in the direction 
perpendicular to the charge current. Thus, when a thin 
FM is deposited on the NM, the pure transverse spin cur- 
rent is injected to the FM and generates the STT. Recent 
experiments [l0l - [T2j demonstrated that the STT gener- 
ated by the spin current can be strong enough to switch 
the magnetization direction of the FM. The magnitude of 
the spin current is often parameterized by the spin Hall 
angle, denoting the ratio between the source charge cur- 
rent density and the resulting spin current density. The 
spin Hall angle was estimated to be 0.06 for Pt and 
—0.15 for /3-Ta. Here, the opposite signs of the Pt 
and /3-Ta spin Hall angles imply opposite spin direction 
of the spin current in the two cases. 

The SHE-induccd STT (SHE-STT) can modify the 
current-driven domain wall (DW) motion. An in-plane 
electric current in a FM/NM system can generate the 
DW motion through two mechanisms. In one mecha- 
nism, the spin-polarized charge current flowing in the 
FM generates adiabatic and nonadiabatic STTs, [Tsl - flTj 
of which effect on the DW motion has been examined ex- 
tensively. The other mechanism is the pure spin current 
injected from the NM due to the SHE. A recent theoreti- 
cal study [111 examined the effect of the SHE-STT on the 
current-driven DW motion in an ideal situation without 
any pinning centers suppressing the DW motion. Such 
an ideal condition is rarely achieved in practical situa- 
tions, however, and it is well known [l9l [2G | that pinning 



centers may affect the DW motion considerably. Also 
from an application point of view, pinning centers are 
needed to achieve reliable control of DW locations in 
DW-based devices. Notches [H or local modulation of 
magnetic properties Q are commonly suggested forms 
of pinning centers for device applications. In this Let- 
ter, we study the SHE-STT effects on the current-driven 
DW motion in the presence of an extrinsic pinning po- 
tential. We show that the SHE can significantly reduce 
the threshold current density to depin a DW from the 
pinning potential, thereby lowering the energy cost for 
the operation of DW-based devices. 

We consider a bi-layered system consisting of a thin 
FM deposited on a NM [Fig. [iKa)]. The current-driven 
DW motion in this system can be described by general- 
izing the standard Landau-Lifshitz-Gilbert (LLG) equa- 
tion d, [H] to include the SHE-STT. In this study, we 
use the collective coordinate approach, (2^, j22l in which 
DW dynamics is described by two collective coordinates, 
DW position q (along cc-direction) and tilting angle (j) 
[see Figs, [ijb) and (c) for its definition in in-plane mag- 
netic anisotropy (IMA) system and perpendicular mag- 
netic anisotropy (PMA) system, respectively]. Although 
this approach has some limitations especially for vortex 
DWs, it was found to be a very useful tool to an- 
alyze transverse DWs that we consider [Figs. [IJb) and 
(c)]. Interestingly, for both IMA [Fig. ^h)] and PMA 
[Fig. [IJc)] systems, this approach results in the identical 
equations of motion describing the coupled dynamics of 
q and 0, [H, [H, [13 
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where a is the Gilbert damping constant, A is the DW 
width, 7 is the gyromagnctic ratio, Hd is the hard-axis 
anisotropy field, bj = [fv^ P / 2eM s) J f is the magnitude 
of the adiabatic STT in the velocity dimension, /36j is the 
magnitude of the nonadiabatic STT, Ms is the saturation 
magnetization, T4xt is an extrinsic pinning potential, and 
Bgjj = TrOsHJN/^tpPJF represents the magnitude of the 
SHE-STT. Here, P and tp are the thickness and spin po- 
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Fig 1: (color online) (a) Schematics of a FM/NM system. 
Magnetic configurations and definition of (j) for a (b) Neel DW 
in IMA system (c) Bloch wall in PMA system. Red arrows 
represent magnetization direction. 

larization of the FM, ^sh is the spin Hall angle of the sys- 
tem, and Jn{Jf) is the current density in the NM(FM). 
Vext is assumed to be a finite-ranged harmonic potential 
placed at (7 = 0: Kxt = - C')0(C " kl) where 

&{x) is the unit step function centered at a; = 0. Here, Vq 
and ^ arc the depth and range of the pinning potential. 
For analysis, it is convenient to rewrite Eq. ([T]) in terms 
of dimcnsionless quantities as 

X-Q(/)=sin2^- J, (2) 

aX + (p^- J(^ + BsHAsin(/)) - VXe{^ - \X\). 

A 

Here, X = q/X, J = 2hjhHd\ V = (2A/A/s-ffdC')K), 
and 6 = {2/-fHd)dO/dt. 

To investigate the SHE-STT effects in the presence of 
V^xt, we first calculate the threshold current density Jc, 
above which a DW can get depinned from the pinning po- 
tential. In our simulation, a DW is placed at the center 
of the pinning potential {X = 0) with the initial tilting 
angle (or polarity) (/>o = or tt. The current pulse is 
then turned on with zero rising time. Figure [2][a) shows 
the numerical simulation result obtained with a = 0.02, 
/3 = 0.01 and £,/X = 1. These parameter choices may be 
applicable to an IMA bi- layer system such as Pt/Py. Q 
In the absence of the SHE [black solid line in Fig.[2l^a)], it 
is well known [l^ that Jc depends on V in three distinct 
ways. In the so-called intermediate regime [0.3 < < 40 
in Fig.dja)], the adiabatic STT SIHJlll is the main de- 
pinning mechanism and Jc is almost independent of V, 
with Jc value similar to the intrinsic threshold value. [l5j 
For smaller (weak pinning regime, V < 0.3) and larger 
(strong pinning regime, V > 40) V, the nonadiabatic 
and adiabatic STTs are respectively the main depinning 
mechanisms overcoming the extrinsic pinning potential 
and Jc increases as V increases. A thick lOnm) 
Py layer 0, H, can have small V in the weak pin- 
ning regime. But when a Py layer is made thin (^ a 
few nm 0), which increases the SHE-STT since Bgu is 
inversely proportional to tp, V is expected to be larger 
due to the enhanced interface contribution to V. Thus 
a SHE bi-layer system is likely to be near the borderline 
between the weak and intermediate pinning regimes. In- 
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Fig 2: (color online) (a) Jc as a function of V for 
BshXcos (fio — —0.6 (blue squares), —0.3 (red circles), 
(black solid line), 0.3 (purple diamonds), and 0.6 (green tri- 
angles). Solid lines connecting symbols are eye-guides. Black 
daslied and dashed-dotted lines represent ( — curves for 
BshXcos (fiQ = —0.3 and —0.6, respectively, (b) X after 
the current injection for Bsh = and J — 0.5 (blue dash- 
dotted line), BsH ~ and J — 0.8 (green dashed line), and 
BshAcos0o = -0.2 and J = 0.5 (red solid line). ( = 0.02 for 
the first two lines and (" — —0.013 for the last line. 



terestingly, it turns out that this is the range where the 
SHE-STT effects are most pronounced. Numerical simu- 
lation of Eq. ^ indicates a significant reduction of Jc in 
this range if i?sHAsgn( J) cos 0o < [—0.3 and —0.6 for 
red circular and blue squared symbols in Fig. [H^a)]. On 
the other hand, Jc changes only mildly in this range if 
SsHAsgn(J) cos(/)o > [0.3 and 0.6 for purple diamond 
and green triangular symbols in Fig.[2l[a)]. Note that for 
a given sign of i?sH (or 6'sh), the SHE reduces Jc sig- 
nificantly only for a DW with a proper sign of cos 0o (or 
proper DW polarity). Hence, one can control the depin- 
ning efficiency by preparing the initial DW polarity to a 
proper value. 

We remark that for metallic ferromagnets with 
IMA our estimation of [BshAI ^ 0.56. Here, we as- 
sume Jn = Jf, ^sh = 0.06; tp = 4nm, P = 0.7, and 
A = 20nm. 

To understand the origin of the Jc reduction, we fur- 
ther analyze (p dynamics. Within the pinning range 
[6(^/A — \X\) = 1], two equations in Eq. ^ can be 
combined into the following 2nd order partial differential 
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Fig 3: (color online) (a), (b), and (c): The DW depinning di- 
rection as a function of J and V for Bsh ~ 0, -BshA cos (jto ~ 
—0.4, and i?sHAcos(/!>o = 0.4, respectively. Orange, gray and 
yellow surfaces represent that the DW is not depinned, the 
DW depins to the left and right, respectively. For the refer- 
ence, electrons flow to the left when J > 0. (d) Jrev (blue 
circles) and Vth (green squares) as a function of _BshAcos(^o. 
Solid and dotted lines connecting symbols are eye-guides. 



equation for (j) dynamics, 

(1 + a^)0+ (2acos2(/) + BshX J cos (f) + aV)(i) (3) 

+ (sin20-J)y = 0. 

For small J, (f> stays near its initial value (f>o, which is 
cither or tt. Thus, to get an insight, we expand above 
equation near <f> = <f>Q to obtain 



'(p -f 2Cwo0 + wg0 = jy/(l + q;2). 



where 



2V 



1 



c 



a(2 + V) + BshA J cose 
v/8F(l + a2) 



(4) 



(5) 



Here, wq and C are an undamped angular frequency and 
the damping ratio, respectively. For B^a = 0, C is always 
positive and 4> approaches towards its stable steady value. 
The accompanied X dynamics is shown in Fig. [2Jb) as 
blue dash-dotted (J < Jc) and green dashed (J > Jc) 
lines. In the presence of the SHE (Bsh 0), however, 
C is not necessarily positive and may become negative if 
SshAJcos(/)o is negative and sufficiently large in magni- 
tude. When C < 0, dynamics does not have a stable 
steady value and results in the DW depinning driven by 
the amplification of (j) [and also X as shown in Fig. [IJb) , 
red solid line]. Thus the condition C = sets an upper 
bound on | Jc|, that is, \Jc\ < Jc,up with 



Jaup = a(2 + F)/|BsHA| 



(6) 



which is shown as black dashed and dash-dotted lines in 
Fig.mja) for _BgHAsgn(J) cos0o = —0.3 and —0.6, respec- 
tively. In the region where the SHE reduces Jc, note that 



Jc,up is in good agreement with numerically calculated 
Jc values [red circles and blue squares in Fig.^a)]. This 
demonstrates that the reduction of Jc by the SHE arises 
because the SHE can induce negative C, and thus amplify 
the (j) dynamics. For i?sHAsgn( J) cos (j)o > 0, on the other 
hand, the SHE increases C and thus Jc as well [purple 
diamond and green triangular symbols in Fig. HJa)]. A 
slight decrease of Jc in the weak pinning regime is due 
to the shift of the stable DW position. Note that one 
sign of BshA J cos (/)o reduces C while the other increases 
it. Thus, to reduce Jc, the DW polarity should be prop- 
erly adjusted before the current injection, for example 
by applying weak in-plane magnetic field or injecting a 
properly designed train of small current pulses. (26j 

Similar reduction of Jc occurs in a PMA bi-layer as 
well, which is natural since IMA and PMA bi-layer sys- 
tems share the same equations of motion for the DW 
depinning dynamics [Eq. ([TJ]. In particular, the upper 
bound Jc,up in Eq. (jB]) applies equally to the PMA bi- 
layer. Thus the difference between the PMA and IMA bi- 
layers may come only through different material parame- 
ters. For the PMA bi-layer system /3-Ta/CoFeB with low 



Gilbert damping a = 0.008 Jc,up is estimated to be 



similar to or even smaller than the corresponding value 
in Fig. ^a) for the IMA bi-layer system. For the PMA 
bi-layer system Pt/CoFeB, the reported value of a ranges 
from 0.025 [ll| to 0.5 [25j. Since Jc,up is proportional to 
a, the reduction of Jc due to the SHE is expected to be 
much less significant for the PMA system with large a. 

Next, we examine the depinning direction. In the ab- 
sence of the SHE (BsH = 0), the DW gets depinned 
to the right (left) of the pinning center when electrons 
flow to the right (left) and thus the depinning direction 
agrees with the electron flow direction [Fig. ^&)] . This 
direction remains the same even in the presence of the 
SHE if BshAJcos^o > 0. On the other hand, when 
SshAJcos^o < and C < 0, interestingly, the depin- 
ning direction is opposite to the electron flow direction 
in a wide J range near Jc [see upper yellow region in 
Fig. El^b) and lower gray region in Fig. EJc)]. 

This situation resembles the situation for the DW mo- 
tion direction in an ideal wire without pinning centers. 
In the absence of the SHE, the terminal velocity of the 
DW is always along the electron flow direction. When the 
SHE exists, however, it was recently demonstrated [l^ 
that the terminal motion direction can be opposite to the 
electron flow direction if Bsk^J cos cf> < 0. Note that for 
both the depinning direction reversal and the terminal 
motion direction reversal, _BshAJcos0 should be nega- 
tive. Thus when the initial DW polarity 0o satisfies this 
condition for a given sign of BshAJ, the DW may get 
depinned against the electron flow direction and main- 
tain the same motion direction. This raises a possibility 
that a DW may propagate against the electron flow di- 
rection without getting pinned in a nonideal nanowire 
containing pinning centers. This may also provide an 
explanation for the reversed DW motion observed ex- 
perimentally, (25I without resorting to the Rashba spin- 
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orbit coupling (RSOC) effect, [l^] For this possibility to 
work, it is iraportant to avoid the DW polarity switch- 
ing or the sign change of cos cj), which sets a constraint; 
J should be smaller than the certain current density 
value (Jrev)j above which the DW polarity switching 
occurs. Thus, J should fall in the current density win- 
dow. Jc < J < ^REV- For such window to exist, the 
inequality Jc < ^rev should be satisfied. By numerical 
simulation, we calculate the threshold pinning potential 
depth Vti^, above which the inequality is violated. Firstly, 
>^REV for V = [blue circular symbols in Fig. ^d)] is 
calculated. For IBshA] = 0.1 ^ 0.6, Jrev is of the or- 
der of 10~^ which amounts to the current density of the 
order o_f lO^^A/m^ for Ms = lO^A/m and Ha = IT. 
Then Jrev is compared with Jc to extract Vth [green 
squared symbols in Fig. |31^d)]. For |-BshA| = 0.1 ^ 0.6, 
Vth is smaller than 0.04 which amounts to 4% of the 
hard-axis anisotropy energy of a DW. [l^] Here, we used 
V = (VoX/Kd^^) and A - ^ where Kd = MsHd/2 is the 
hard-axis anisotropy energy. To observe the DW motion 
against the electron flow direction, V < Vth should be 
satisfied. Thus to explain the reversed DW motion in a 
recent experiment 12511 with the SHE-STT (but without 
the RSOC effect [23|), pinning centers should be very 
weak. 



Lastly, recent theory [131 indicates that the RSOC- 
induced STT (RSOC-STT) may drive the DW against 
the electron flow direction in a FM sandwiched by two 
dissimilar layers such as metal and oxide. In this case, 
the so-calle field-like STT arising from the RSOC sup- 
presses the DW polarity switching. Thus the DW motion 
against the electron flow direction can be maintained for 
relatively larger value of V compared to the case with 
the SHE-only (without RSOC-STT) case. The strength 
of the extrinsic pinning potential can be controlled, for 
example by introducing the notch in the wire. Thus this 
feature may be used to distinguish the RSOC-STT and 
the SHE-STT in experiments. 

In summary, we examine the SHE-STT effects on a 
DW motion in the presence of the extrinsic pinning po- 
tential. The presence of the SHE-STT mainly affects the 
damping ratio of a DW precession in the pinning po- 
tential and by inducing negative damping ratio, it can 
signiflcantly reduce the threshold current density. We 
also examined the DW motion direction and found that 
the SHE-STT can induce a DW motion against the elec- 
tron flow direction if the pinning potential is sufficiently 
weak. By using this dependence on the pinning potential 
strength, we suggest one way to distinguish the SHE-STT 
and RSOC-STT in the presence of the pinning potential. 
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